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Abstract 

We consider the evolution of a parton plasma created in Au+Au collisions at 
LHC and at RHIC energies. Using Boltzmann equation, relaxation time approxi¬ 
mation and perturbative QCD, we show the physics of both thermal and chemical 
equilibration in a transparent manner. In particular, we show inelastic processes are, 
contrary to common assumption, more important than elastic processes, the state 
of equilibration of the system can compensate for some powers of a s for the pur¬ 
pose of equilibration, the two-stage equilibration scenario is, barring any unknown 
non-perturbative effects, inevitable and is an intrinsic feature of perturbative QCD, 
and gluon multiplication is the leading process for entropy generation. 


1 Introduction 

Essential questions to ask in heavy ion collisions, especially at RHIC and at LHC, as¬ 
suming a gas of weakly interacting partons can be formed in the central region, is how 
fast will this quark and gluon system approaches equilibrium if the expansion is not too 
rapid for the interactions. If indeed this parton gas can approach equilibrium, will it be 
able to end up as a quark-gluon plasma? And what is the degree of equilibration can one 
reasonably expect when the phase transition sets in? In short, is equilibration fast and 
how fast? 

A number of previous works have already attempted at providing answers to these and 
other related questions of equilibration. However, due to the fact that non-equilibrium 
problems are difficult, they addressed either only thermalization SI,SB or only parton 
chemical equilibration || [|]. Furthermore, previous attempts at the thermalization prob¬ 
lem have only been done in a heuristic manner without really using QCD interactions. 
With the exception of the parton cascade model (PCM) [|7j, which is based on present 
knowledge of perturbative QCD and some very involved computations, only then it is able 
to consider both thermalization and chemical equilibration simultaneously as it should be 
and as it happens in the collisions. In this talk, we present a relatively simple way to do 
this and hence keeping things simple and clear so that the physics becomes transparent. 

A much used assumption and starting point in the studies of the physics of heavy ion 
collisions is kinetic equilibration is rapid <, 1.0 fm/c and hydrodynamics expansion is 
well underway. PCM has shown that such short rapid thermalization is too optimistic. In 
the following, we will consider the evolution of a parton plasma and show that complete 
kinetic and chemical equilibration are slow, although interactions can indeed dominate 
over the expansion, and full equilibration cannot be attained before the phase transition 
assumed to be at T c ~ 200 MeV. We also show that inelastic processes are more important 
than elastic ones in equilibration contradicting the common untested assumption of the 
contrary. 
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2 To Determine the Evolution of a QCD Plasma 


To study the evolution of a plasma of quarks and gluons, it is sufficient to know the 
particle distributions. For this purpose, we use the set of rather standard assumptions for 
relativistic heavy ion collisions and Baym’s form of the Boltzmann equation |l| 
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where r = \/t 2 — z 2 , with the collision terms on the right hand side approximated by the 
relaxation time approximation 
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where f eq is the equilibrium distribution f eq = l/(exp(p/T eq ) =F 1) and 9{r) is the collision 
time. Whether the plasma equilibrates or not and how fast does it equilibrate depends 
very much on 6 { t ) [f|, ||. 

With the combination of Eq. (jl]) and Eq. one can already write down a solution 
to Eq. ([I]), which depends, however, on the two numerical parameters T eq and 9 that need 
to be determined from QCD. We use the simplest QCD interactions at the tree level for 
the collision terms for this purpose. They are 
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Here to keep things simple and for our purpose, it is sufficient to include the simplest two 
leading inelastic processes^ the first term of Eqs. (||) and (|4|) and all the binary elastic 
processes. Quarks and gluons will be treated as different particle species and not as 
generic partons so that their respective distributions, f g and f q , are governed by different 
Boltzmann equations as they should be and depend on different collision times, 9 g and 
9 q , and different equilibrium temperatures, T eqg and T eqq . 

Using perturbative QCD and suitably infrared regularized the matrix elements by 
medium effects, one can construct the collision terms for gluons C g and for quarks C q 
semi-classically in the usual way [Q. For soft gluon emissions, Landau-Pomeranchuk- 
Migdal effect |J has to be incorporated due to multiple scatterings in the medium || || . 

With the real collision terms from QCD, one can construct two equations for each 
particle species to solve for the two time-dependent unknowns in the particle distribution. 
We choose the following rate equations. 

1) Energy density rate in an one-dimensional expanding system 
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2 Here we assign a single chemical potential, n q , to all rif flavours of fermions so flavour changing 

interactions in the first interactions of Eq. (§ is not considered as inelastic. 
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2) Collision entropy density rate 
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where i = g,q, q, Pl is the longitudinal pressure defined later in Eq. (||) and u t is the 
multiplicity, for gluons n g = 2x8 and for quarks n q — 2 x 3 x n/. These equations are 
constructed from Eqs. © and (||) and Cf s are now understood to be the real QCD collision 
terms. From them, dfs and T eq f s can be determined. The collision entropy equations, 
in fact, allow one to break down the entropy generation process due to collisions into 
each of its contributing elements and find out which processes are more important for 
equilibration. We will show this later on in Sect. 


3 The Approach to Equilibrium 


We take the initial conditions at the isotropic moment To from HIJING |T0| results for 
Au+Au collisions at RHIC and at LHC The various equations described in Sect. 

|H are solved numerically and the various collective variables can be calculated. The 
numerical details and parameters for solving the equations can be found in ||. Here we 
concentrate on the results and the physics of equilibration. 

As mentioned in the introduction, there are two types of equilibration: chemical 
and thermal, which happens simultaneously in heavy ion collisions. To check the par- 
ton composition in the plasma, it is common to use the concept of fugacity defined 
as f = exp (fa/T) so that the kinetically equilibrated distributions can be written as 
fi = l/(exp(p/T i )/h 1 =F !)• Since chemical equilibration is essentially fa ^ 0 starting from 
some initial negative values so that in terms of fugacity, /, approaches 1.0. Clearly Ifs 
as well as T's exist only when local kinetic equilibrium has been achieved which is not 
the case most of the time during the evolution of the parton gas. Nevertheless, one can 
estimate these quantities from the energy densities e/s and number densities nf s in the 
usual way [||. 

These are shown in Fig. |l]. The solid lines are for gluons and the dashed lines are 
for quarks. The evolution is stopped when the temperature estimates all fall below 200 
MeV, the assumed phase transition temperature. As can be seen, the gluon fugacities 
approach 1.0 much more rapidly than those of the quarks both at RHIC and at LHC. 
Gluon fugacities are close to 1.0 near the end as a result but not those of the quarks as 
in agreement with previous studies [| I- §■ 

To check for kinetic equilibration, there is not the equivalent quantity of as for 
chemical equilibration so instead one checks the isotropy of momentum distribution by 
comparing the ratio of the longitudinal to transverse pressure Pl/pt and the ratio of a 
third of the energy density e/3 pp to the transverse pressure. These pressures are defined 
by 

d 3 ~ 


i \ [dp P z x , >, 

PL,T i{ T ) = Vi / fiiP±,Pz,T) 

J (2n r v 
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These quantities are related of course by e = 2 pp + Pl to the energy density, so isotropy 
means pr = Pl = e/3. These ratios should approach 1.0 in an equilibrating parton gas. 
We have plotted these ratios in Fig. § for gluons and for quarks. 
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Figure 1: The evolution of the species temperatures T g and T q and the fugacities l g and 
l q at LHC and at RHIC. The solid (dashed) lines are for gluons (quarks). The rising 
(falling) curves are the fugacity (temperature) estimates. 
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Figure 2: The plots of the evolution of the ratios of the longitudinal pressure (solid line) 
and one-third of the energy density (dashed line) to the transverse pressure, Pl/Pt and 
e/3 pt respectively, for gluons and for quarks at LHC and at RHIC. 


The ratios are 1.0 initially because we have started from an isotropic momentum 

configuration. As the system expands, the distributions become anisotropic and reach 

maximum anisotropy quickly. The subsequent return towards isotropy is only progressive 
and at the end, complete isotropy is not fully recovered. For the gluons, one can perhaps 
argue an approximate isotropy has been achieved near the end especially at LHC. For 
the quarks, they have not yet passed the half-way mark. To check that, although slow, 
these are indeed thermalization behaviours, we can compare with the case that the sys¬ 
tem ends up in free streaming. In that case, interactions are not important and can be 
simulated by letting 6 g and 6 q to oo and the particle distributions are described by / 0 . 
The corresponding pressures and ratios work out, as r —> oo, to be 

Pl/Pt —> 2tq/t 2 — > 0 and e/3p T — > 2/3 (9) 

which are clearly not the behaviours in Fig. ^j. This provides clear evidence that the 
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Figure 3: The evolution of the collision entropy density ratios of gg <—> gg (dashed), 

gg < -> qq (solid), and gq < -> gq or gq <-> gq (dot-dashed) to gg < -> ggg for gluons 

(top figures) and gq < -> gq or gq < -> gq (solid) and the sum of all fermion clastic 

scatterings to gg —> qq for quarks (bottom figures). 

interactions indeed dominate over the expansion. 

Finally, we show the dominant processes in the equilibration of the plasma. It is 
common to assume that thermalization is driven by elastic processes in the studies of the 
various physics in heavy ion collisions. Inelastic processes are relegated to the minor role 
of essentially only for chemical equilibration. To show that this is not true, we break 
down the gluon and quark collision entropy density rate, ( ds g /dr) co u and ( ds q /dT) co u 
respectively, to their contributing elements and plot their ratios. The ratios plotted in 
Fig. ^ are each process to gluon multiplication for gluon and to gluon-gluon conversion 
into quark-antiquark pair for quark. 

Initially, all the ratios are below 1.0 in the top figures i.e. gluon multiplication is 

dominant, at some point around 2.0 fm/c at LHC and 4.0 fm/c at RHIC. gg <-> qq 

rises above one (solid line) and overtakes gg <—> ggg as the dominant process. The 
elastic scattering ratios (dashed and dot-dashed lines) remain below 1.0. In the bottom 

figures, those for quarks, gg < -> qq remains dominant throughout. As can be seen, all the 

curves remain below 1.0. So for the gluons, gluon multiplication dominates initially, when 
this starts to slow down as equilibrium is getting near and gluons are near saturation, 
gluon-gluon annihilation into quark-antiquark pair and the reverse process become more 
important because quarks and antiquarks are still far from completing the equilibration. 
For this reason, in the top figures, gluon-fermion elastic scattering ratios continue to rise 
but that of the gluon-gluon elastic scattering is not changing very much with r. So it 
is clear that inelastic processes are more important contrary to common assumption and 
that gluon multiplication leads in the production of entropy as long as they are still far 
enough from full equilibrium. 

A reason that intuitively elastic processes should be more important is because of 
the larger cross-sections. Inelastic processes such as the radiation of an extra gluon off 
a two-body elastic scattering is down by a s for instance. This reasoning is, however, 
incorrect. In a medium, any interaction that can happen can also happen in the reverse 
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direction, it is not simply the forward or backward reaction that enters the collision terms 
of the Boltzmann equation Eq. ({!]) but the difference of the two. So it is not the sizes of 
the interaction cross-sections which determine what processes should be or should not be 
more important in the equilibration of a many body system. 

To summarize, in this talk, we have shown that strictly speaking equilibration cannot 
be completed in heavy ion collisions at RHIC and at LHC. One can at best consider the 
system as a fluid mixture of an approximately therm a h zed and chemical equilibrated gluon 
plasma and a still far from equilibrated quark and antiquark plasma. The usual belief that 
elastic processes are responsible for thermalization is flawed. As we have shown, inelastic 
processes are even more important. One has to be very careful in comparing interactions 
in a medium, if it is done by relying on the sizes of the scattering cross-sections alone, 
one is prone to getting the wrong answer. 
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